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We present measurements of the B — > nK-y branching fractions and upper limits for the B — > rj'K'y 
branching fractions. For B + — + r]K +r y we also measure the time-integrated charge asymmetry. 
The data sample, collected with the BABAR detector at the Stanford Linear Accelerator Center, 
represents 232 x 10 6 produced BB pairs. The results for branching fractions and upper limits at 
90% C.L. in units of 10~ 6 are: B(B° -> V K ~f) = 11.3±|;| ± 0.6, B(B+ -» 77^+7) = 10.0 ± 1.3 ± 0.5, 
B(B° -> rj'K°-y) < 6.6, B{B + -> rj'K+'y) < 4.2. The charge asymmetry in the decay B + -» 77^+7 
is A c h = —0.09 ± 0.12 ± 0.01. The first errors are statistical and the second systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Radiative B meson decays have long been recognized 
as a sensitive probe to test the Standard Model (SM) 
and to look for new physics (NP) [1| . In the SM, flavor- 
changing neutral current processes such as b — > S7 pro- 
ceed via radiative loop (penguin) diagrams. The loop dia- 
grams may also contain new heavy particles, and there- 
fore are sensitive to NP. Measurements of the branch- 
ing fractions of a few of the exclusive decay modes ex- 
ist: if*l892) 7 Hi, Ki(1270>y 0, #$(1430)7 [10, 
Kmr^f 5], <\>K^ [6j and K + r] r ) |7(. The measured branch- 
ing fraction of inclusive b — > sj and exclusive radiative B 
decays are in good agreement with SM predictions . 
Direct ancl mixing-induced [ll| CP asymmetries in 
exclusive radiative B decays are expected to be very 
small in the SM. Measurement of direct CP asymmetries 
in exclusive radiative decays can provide a clear sign of 
NP [13. 

We present analyses of the exclusive decay modes 
B + — * r/K + j and B° — > r/K°j 13] , which have pre- 
viously been measured by the Belle Collaboration 0, 
and B + — > rj'K + j and B° — > i]'K°j which are studied 
for the first time. The results presented here are based 
on data collected with the BABAR detector a t the 
PEP-II asymmetric-energy e + e~ collider 0] located at 
the Stanford Linear Accelerator Center. The analyses 
use an integrated luminosity of 211 fb -1 , corresponding 
to 232 million BB pairs, recorded at the Y(4S) resonance 
(at a center-of-mass energy of ^/s = 10.58 GeV). 

Charged particles from e + e~ interactions are detected, 
and their momenta measured, by a combination of a ver- 
tex tracker (SVT) consisting of five layers of double-sided 
silicon microstrip detectors, and a 40-layer central drift 
chamber (DCH), both operating in the 1.5 T magnetic 
field of a superconducting solenoid. We identify photons 
and electrons using a CsI(Tl) electromagnetic calorime- 
ter (EMC). Further charged-particle identification is pro- 
vided by the average energy loss {dE/dx) in the track- 
ing devices and by an internally reflecting ring-imaging 
Cherenkov detector (DIRC) covering the central region. 
A K/tt separation of better than four standard deviations 
is achieved for momenta below 3 GeV/c, decreasing to 
2.5(7 at the highest momenta in the B decay final states. 



We reconstruct the primary photon, originating from 
the B decay candidate, using an EMC shower not associ- 
ated with a track. We require that the photon candidate 
fall within the fiducial region of the EMC, has the ex- 
pected lateral shower shape, and is well-separated from 
other tracks and showers in the EMC. The primary pho- 
ton energy, calculated in the T(45) frame, is required to 
be in the range 1.6 - 2.7 GeV. We veto photons from 
TT°(r]) decays by requiring that the invariant mass of the 
primary photon candidates combined with any other pho- 
ton candidate of laboratory energy greater than 50 (250) 
MeV not be within the range 115-155 (507-587) MeV/c 2 . 
Charged K candidates are selected from tracks, by us- 
ing particle identification from the DIRC and the dE / dx 
measured in the SVT and DCH. 

The B decay daughter candidates are reconstructed 
through their decays 7r° — > 77, 77 — > 77 (r? 77 ), 77 — > 
TT+ir-ir (773^), 77' -> ?7 77 7r+7r- (r}' v „ n ), and 77' -> p°j 
(r]'), where p° — > 7T+7T - . Here we require the laboratory 
energy of the photons to be greater than 50 MeV (200 
MeV for 7]'). We impose the following requirements 
on the invariant mass in MeV/c 2 of these particles' final 
states: 120 < 771(77) < 150 for vr°, 490 < 771(77) < 600 
for 77 77 , 520 < m(ir + iT-iT Q ) < 570 for 773^, 930 < 
m{iT + TT-ri) < 990 for 77^, 910 < 777(^+^-7) < 1000 
for 77^, and 510 < m(7r+vr-) < 1000 for p°. For the 
77' and 77 these requirements are sufficiently loose as to 
include sidebands, since these observables are used in 
the maximum-likelihood (ML) fit described below. Sec- 
ondary pions in 77' and 77 candidates are rejected if their 
DIRC and dE/dx signatures satisfy tight requirements 
for being consistent with protons, kaons, or electrons. 

Neutral K candidates are formed from pairs of 
oppositely-charged tracks with a vertex \ 2 probability 
larger than 0.001, 486 < m(7r+7r") < 510 MeV/c 2 and a 
reconstructed decay length greater than three times its 
uncertainty. We require the momentum of the 77 or 77' 
in the T(4S) frame to be greater than 0.9 GeV/c (0.6 
GeV/c in modes with 77' „.). The invariant mass of 77 AT 
and 77'^ systems is required to be less than 3.25 GeV/c 2 . 
In rj'K^/ final states, we suppress background from the 
decay J/tpK, with J/tp — * 77 ; 7 by applying a veto on the 



5 



reconstructed rj'j invariant mass. Denning the helicity 
frame for a meson as its rest frame with polar axis along 
the direction of the boost from the parent rest frame, and 
the decay angle #dcc as the polar angle of a daughter mo- 
mentum in this helicity frame, we require for the rf de- 
cays I cos9 p dcc \ < 0.9, and for ?7 77 decays | cosfi^J < 0.9, 
to suppress combinatorial background. 

A B meson candidate is reconstructed by combin- 
ing an rj or rf candidate, a charged or neutral kaon 
and a primary photon candidate. It is characterized 
kinematically by the energy-substituted mass toes = 
\J (s/2 + p • Pb) 2 /-Eq — P% an d energy difference AE = 
E* B — \yfs, where the subscripts and B refer to the 
initial T(45 l ) and to the B candidate in the lab-frame, 
respectively, and the asterisk denotes the T(45) frame. 

Background arises primarily from random track com- 
binations in e + e~ — * qq events. We reduce this back- 
ground by using the angle 9t between the thrust axis 
of the B candidate in the T(4S I ) frame and the thrust 
axis of the rest of the event. The distribution of | cos#t| 
is sharply peaked near 1 for combinations drawn from 
jet-like qq events, and is nearly uniform for BB events. 
We require |cos$t| < 0.9. Furthermore events should 
contain at least the number of charged tracks in the can- 
didate decay mode plus one. For ri 11 K + ^i we require 
at least 3 charged tracks in the event. If an event has 
multiple B candidates, we select the candidate with the 
highest B vertex \ 2 probability. 

We obtain signal event yields from unbinncd extended 
maximum-likelihood fits. The principal input observ- 
ables are AE, toes and a Fisher discriminant T . Where 
relevant, the invariant masses TO res of the intermediate 
resonances and | cos ? | are also used. The Fisher dis- 
criminant T combines four variables: the angles with 
respect to the beam axis of the B momentum and the 
B decay product's thrust axis (in the T(45) frame), and 
the zeroth and second angular moments Lo,2 of the en- 
ergy flow about the B° thrust axis. The moments are 
defined by Lj — J2iPi x |cos#i| J , where 9t is the angle 
with respect to the B thrust axis of track or neutral clus- 
ter i, pi is its momentum, and the sum excludes the B 
candidate daughters. 

We estimate BB backgrounds using simulated samples 
of B decays [lij . Signal and inclusive b — > events are 
simulated according to the Kagan-Neubert model [l7| . 
The branching fractions in the simulation are based on 
measured values and theoretical predictions. 

For each event i and hypothesis j (signal, continuum 
or BB background), the likelihood function is 



N 



£ = e -(E »,) J] 



3 

E 

3=1 



(1) 



where N is the number of input events, rij is the num- 
ber of events for hypothesis j and Vjfa) is the corre- 
sponding probability density function (PDF), evaluated 



with the observables of the ith event. Since corre- 
lations among the observables are small, we take each 
V as the product of the PDFs for the separate vari- 
ables. We determine the PDF parameters from Monte 
Carlo simulation for the signal and BB background, 
while using sideband data (5.25 < toes < 5.27 GeV/c 2 ; 
0.1 < \AE\ < 0.2 GeV) to model the PDFs of contin- 
uum background. We parameterize each of the func- 
tions "Psig(mEs), VsigiAE), Vj{T), and the compo- 
nents of Vj(m rcB ) that peak in toes with either a Gaus- 
sian, the sum of two Gaussian distributions, or an asym- 
metric Gaussian function, as required, to describe the 
distribution. Distributions of AE for BB and contin- 
uum background and | cos #j ec | are represented by linear 
or quadratic functions. The BB and continuum back- 
ground in toes is described by the ARGUS function 
xyl — x 2 ex p [— £(1 — x 2 )] , with x = 2toes/\A and a 
parameter £ [l8j. We allow continuum background PDF 
parameters to vary in the fit. 

Large control samples of B decays to charmed final 
states of similar topology are used to verify the simu- 
lated resolutions in AE and toes- Where the control data 
samples reveal differences from the Monte Carlo (MC) in 
mass or energy resolution, we shift or scale the resolu- 
tion used in the likelihood fits. Any bias in the fit is 
determined from a large set of simulated experiments in 
which the qq background is generated from the PDFs, 
and into which we have embedded the expected number 
of BB background and signal events chosen randomly 
from fully simulated Monte Carlo samples. 

In Table [I] we show the signal yield, the efficiency, and 
the product of daughter branching fractions for each de- 
cay mode. The efficiency is calculated as the ratio of the 
number of signal MC events entering into the ML fit to 
the total generated. We compute the branching fractions 
from the fitted signal event yields, reconstruction effi- 
ciencies, daughter branching fractions, and the number 
of produced B mesons, assuming equal production rates 
of charged and neutral B pairs. We correct the yield 
for any bias estimated by the simulations. We combine 
results from different channels by combining their like- 
lihood functions, taking into account the correlated and 
uncorrelated systematic errors. We report the statistical 
significance and branching fraction for the individual de- 
cay channel; for combined measurements having a signifi- 
cance smaller than 5cr, we also report the 90% confidence 
level (C.L.) upper limit. 

The statistical error on the signal yield is taken as the 
change in the central value when the quantity —2 InC 
increases by one unit from its minimum value. The sig- 
nificance is the square root of the difference between the 
value of —2 In £ (with systematic uncertainties included) 
for zero signal and the value at its minimum. The 90% 
C.L. upper limit is taken to be the branching fraction 
below which lies 90% of the total likelihood integral in 
the positive branching fraction region. 
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TABLE I: Signal yield, detection efficiency e (%), daughter branching fraction product 11^*1 significance S (a) (including 
systematic uncertainties), measured branching fraction B with statistical error for each decay mode. For the combined mea- 
surements we give the significance (with systematic uncertainties included) and the branching fraction with statistical and 
systematic uncertainty (in parentheses the 90% C.L. upper limit). For the r]K + 'y mode we also list the measured signal charge 
asymmetry A c h- 
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The measured charge asymmetry in the decay 
B + — > r]K + j is corrected for an estimated bias of 
—0.005 ± 0.010, determined from studies of signal Monte 
Carlo events and data control samples and from calcula- 
tion of the asymmetry due to particles interacting in the 
detector. The result is A ch = -0.09 ± 0.12 ± 0.01 with 
an asymmetry interval [-0.282,0.113] at 90% C.L.. 

Figure Q] shows, as representative fits, the projections 
onto toes and AE for the decays r]K + 'y, rjK°j, rj'K + j 
and rj'K°j for a subset of the data for which the signal 
likelihood (computed without using the variable plotted) 
exceeds a mode-dependent threshold that optimizes the 
sensitivity. 

Figure shows the distribution of the r]K invariant 
mass for signal events obtained by the event-weighting 
technique (sPlot) described in Ref. 01 . We use the co- 
variance matrix and PDFs from the ML fit to determine 
a probability for each signal event. The resulting distri- 
butions (points with errors) are normalized to the signal 
yield. This mass distribution is useful to compare with 
theoretical predictions for radiative decays [l7|. 

The main sources of systematic error include uncer- 
tainties in the PDF parameterization and ML fit bias. 
For the signal, the uncertainties in PDF parameters are 
estimated by comparing MC and data in control sam- 
ples. Varying the signal PDF parameters within these 
errors, we estimate yield uncertainties of 1-2 events, de- 
pending on the mode. The uncertainty from fit bias is 
taken as half the correction itself (1-3 events). System- 
atic uncertainties due to lack of knowledge of the pri- 
mary photon spectrum are estimated to be in the range 
2-6% depending on the decay mode. Uncertainties in our 
knowledge of the efficiency, found from auxiliary studies, 
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N t and 1.5% x N~, where N t and 



FIG. 1: The B candidate m^s and AE projections for r/K + 'y 
(a, b), r)K°-y (c, d), Tj'K+'y (e, f) and rj'K ^ (g, h). Points 
with error bars represent the data, the solid line the full fit 
function, and the dashed line its background component. 



the B candidate. There is a systematic error of 2.1% 
in the efficiency of K® reconstruction. The uncertainty 
in the total number of BB pairs in the data sample is 
1.1%. Published data |8j provide the uncertainties in the 
B daughter product branching fractions (0.7-3.4%). We 
assign a systematic uncertainty of 0.010 to A c h for the 
bias correction. 

In conclusion, we have measured the central values 
and 90% C.L. upper limits in units of 10~ 6 for the 



are the numbers of tracks and photons, respectively, in branching fractions: B(B° — > rjK°"f) = 11.3T 2 'g ± 0.6 
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FIG. 2: Plot otrjK invariant mass for signal using the weight- 
ing technique described in the text for the combined sub-decay 
modes: (a) charged modes, (b) neutral modes. 

B(B+ -> r)K+i) = 10.0 ± 1.3 ± 0.5 , B{B° -> rj'K°j) = 

l.itti ± o- 1 (< 6 - 6 )> B ( B+ -> v'k + i) = i-9±i:i ± 

0.1 (< 4.2). The measured branching fractions of the 
decay modes B + — > r]K + ^ and i3° — > r)K ^ are in good 
agreement with the values reported by the Belle Col- 
laboration We do not find evidence of the decays 
B" -> 77'^% and B+ r/Jf+7. The S -> decays 
may be suppressed with respect to B — > 77X7 decays due 
to destructive interference between two penguin ampli- 
tudes [2(J. This effect has been measured in B decays to 
T]' K and rjK, for which the branching fraction of the for- 
mer is enhanced with respect to that of the latter |2l| . We 
have also measured the charge asymmetry in the decay 
B+ -► rjK+j to be A ch = -0.09 ± 0.12 ± 0.01, consistent 
with zero. The A ch interval at 90% C.L. is [-0.28,0.11]. 
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